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For an extended Harper model, the fidelity for two lowest band edge states corresponding to dif- 
ferent model parameters, the fidelity susceptibility and the von Neumann entropy of the lowest band 
edge states, and the spectrum-averaged von Neumann entropy are studied numerically, respectively. 
The fidelity is near one when parameters are in the same phase or same phase boundary; otherwise 
it is close to zero. There are drastic changes in fidelity when one parameter is at phase boundaries. 
For fidelity susceptibility the finite scaling analysis performed, the critical exponents a, /3, and v 
depend on system sizes for the metal-metal phase transition, while not for the metal-insulator phase 
transition. For both phase transitions vja, ~ 2. The von Neumann entropy is near one for the 
metallic phase, while small for the insulating phase. There are sharp changes in von Neumann 
entropy at phase boundaries. According to the variation of the fidelity, fidelity susceptibility, and 
von Neumann entropy with model parameters, the phase diagram, which including two metallic 
phases and one insulating phase separated by three critical lines with one bicritical point, can be 
completely characterized, respectively. These numerical results indicate that the three quantities 
are suited for revealing all the critical phenomena in the model. 

PACS numbers: 71.30.+h, 03.67.Ud, 71.23.Ft 



I. INTRODUCTION 



In recent years, tools from the quantum-information 
theory 1 '-, specifically the ground state fidelity 3 and quan- 
tum entanglement^, have been widely exploited to 
characterize quantum phase transitions(QPTs) 6 . For 
example, in one dimensional XY and Dicke model, 
fidelity between two ground states corresponding to 
slightly different values of the parameters drastically de- 
creases at phase transition points 3 . Subsequently, sim- 
ilar properties are also found in fermionio 7 -^, bosonic 
systems 9 -^ and other various spin system s 11 ' 12 . Very 
recently, fidelity susceptibility(FS) (the second deriva- 
tive of fidelity) is introduced to signal QPTs in 
one-dimensional Hubbard model s 13 ' 14 ! 15 , the Lipkin- 
Meshkov-Glick models, the Kitaev honeycomb modeli 7 -, 
and various spin system s 18 ' 19 ' 20 . It is found that FS is 
more curial than fidelity itself for it does not depend 
on the slightly difference values of model parameters. 
In Refsi 21 ' 22 , the fidelity between arbitrary two ground 
states is studied in one dimensional quantum Ising model. 
Singularities are found in fidelity surfaces for QPTs 22 . 
The main advantage of the fidelity to identify QPTs is 
thati^, it need not a priori knowledge of the order param- 
eter, topology, etc, since the fidelity is a purely Hilbert- 
space geometrical quantity. 



At the same time, quantum entanglement has 
been extensively applied in condensed matter 
physics. 23 ' 24 ' 25 ' 26 ' 27 ' 28 ' 29 For example, quantum en- 
tanglement measured by the von Neumann entropy 
has been studied in the Hubbard model for the dimer 
case 2 ^, in the extended Hubbard model for different 



band filling 2 ^, in quantum small-world networks 2 ^, and 
in low-dimensional semiconductor systems^ 6 -. It is found 
that the von Neumann entropy is suitable for analyzing 
the interplay between itinerant and localized features^, 
as well as characterizing quantum phase transitio n 24 ' 27 
and the localization-delocalization transition of electron 
states^s^. 

On the other hand, since the Hofstadter butterfly en- 
ergy spectrum was found in 19762°-, the problem of elec- 
trons in two-dimensional periodic potential in a magnetic 
field has been attracted much attention^^&^ii^ 5 .^ . 
After fixing the quasimomentum in one of the directions, 
a one-dimensional quasiperiodic system called the Harper 
model is deduced^ .. The system shows interesting metal- 
insulator transitions(MIT) 32,33 . Considering the next- 
nearest-neighbor hopping for electron on the square lat- 
tice in a uniform magnetic field, an extended Harper 
model are proposed^ 4 - and studied extensively^ 5 .. Very 
recently, a similar extended Harper model is introduced 
from two-dimensional electrons on the triangular lattice 
in a uniform magnetic field 36 . The phase diagram has a 
very rich structure, which shown in FigQ] In regions I 
and III the wave functions (spectra) are extended (abso- 
lutely continuous), and in the region II the wave func- 
tions are localized (pure points). On the three boundary 
lines, the wave functions (spectra) are critical (singular 
continuous). Besides the traditional MIT, there are novel 
transitions between the two metallic phases(MMT). At 
the bicritical point where the triangular lattice symmetry 
is retained, both level statistics and multifractal analysis 
show quantitively different behaviors from those of other 
critical points^ 6 -. 

Considering the above two aspects, we perform a de- 
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FIG. 1: Phase diagram of the extended Harper model. 



tailed study of the fidelity between arbitrary two quan- 
tum states, FS and von Neumann entropy for the ex- 
tended Harper model^ 6 -. For each of the three quanti- 
ties, there are drastic changes at phase boundaries, i.e., 
the phase diagram can be distinguished according to the 
variations of them with model parameters. Our studies 
provide that the two tools, fidelity and von Neumann en- 
tropy, borrowed from the quantum-information theory, 
are well enough to identify phase transitions in the sys- 
tem. 

The paper is organized as follows. In the next section 
the extended Harper model and the definitions of fidelity, 
FS and von Neumann entropy are introduced. In Sec. lIIII 
the numerical results are presented. And we present our 
conclusions and discussions in Section IIVI 



of an electron with eigenenergy E 1 is 

i*7> = E^n = E^ c «i°)' ( 2 ) 

n n 

where ip^ is the amplitude of the 7th eigenstate at the 
nth site. If set A — 2^, /1 — I s - and t a is taken as units, 
the eigenvalue equation 36 becomes 

-[1 + ^e _27r ^(" _1 / 2 ' +ifey ]T/' n -i — [1 + ^te 27 ™ < ^ n+1 / 2 )~ lfcM ] 
ip n+ i - A cos(27T0n + k y )tp„ = Eip n . (3) 

At fi = and <p is irrational, this is reduced to the 
Harper equation. Intensively analytical and numerical 
studie s 32 ' 33 for the Harper model show that for A > 2 
the spectrum is pure-point like and all eigenstates are ex- 
ponentially localized. For A < 2 the spectrum becomes 
continues with dclocalizcd eigenstates corresponding to 
ballistic classical motion. For A = 2 the situation is criti- 
cal with a singular-continuous multifractal spectrum and 
power law localized eigenstates. MIT can occur at A = 2. 

B. Fidelity 

Let l^/oCA,/^)) denote the lowest band edge state. Ac- 
cording to the definition of fidelity-^ ' 21 ' 22 , the quantum 
fidelity (or the modulus of the overlap of eigenstates) is 
given by 

F(A,/x;A ,Mo) = |(*o(A,/x)|*o(A ,Mo)>|. (4) 
Obviously, F — 1 if A = Aq and [i = /io- 



C. Fidelity susceptibility 



II. THE EXTENDED HARPER MODEL, 
FIDELITY, FIDELITY SUSCEPTIBILITY AND 
VON NEUMANN ENTROPY 

A. The extended Harper model 

The tight-binding Hamiltonian for an electron moving 
on a triangular lattice in a magnetic field^ 6 - can be re- 
duced to 

H= - J2l t a+ t c e ~ 2ni ' l ' (n ~ 1/2)+iky ]<i C n-l 
n 

n 

- 2^4 cos(27n/)n + fcj,)c^c„, (1) 

n 

where t a , t& and t c are the hopping integral for each bond 
on the triangular lattice, (/) is a flux that a uniform mag- 
netic field penetrates each triangle, k y is a momentum 
in the y direction, a n (c„) is the creation (annihilation) 
operator of the nth site in the x direction. 

Let \n) denote |0, 0), the general eigenstate 



Similarly as that shown in Ref^, the fidelity for two 
lowest edge states with a slightly different parameter val- 
ues is defined as 



F{q) = \^ a (q)\^ a {q + 5q))\. 



(5) 



For simplicity, a certain path q = g(A, /x) in parameter 
spaces can always be supposed. Then the FS can be 
calculated a a 13 ' 17 ^ 



Xf = lim 



5q->0 Sq 



2\n F(q) ^ „ h 

-= E 9a, b n a n b , (6) 



a—\,fi;b—X,^i 



where n A = dq/dX (n M = dq/dfi) denotes the tangent 
units vector at the give parameter point (A,/x). For the 
present model, let define the driving Hamiltonians 

H\ = — cos(27T(/>n + kyjc^Cn (7) 



and 



n 

J2[e 2 " Hn+1/2) - lky ]clc n+1 . (8) 
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D. von Neumann entropy 



(9) 



The general definition of entanglement is based on the 
von Neumann entropy 3 -^. For an electron in the system, 
there are two local states at each site, i.e., |0) n , |l) n . The 
local density matrix p n is define d 23 ' 24 ! 25 ' 29 by 

Pn = z n |1> „„ (1| + (1 - z n ) |0) „„ (0| , (10) 

where z n — (^ 7 | c^Cn |\t 7 ) = \ipZ\ 2 1S the local occupation 
number at the nth site. Consequently, the corresponding 
von Neumann entropy related to the nth site is 



E 1 



-z n \og 2 z n -(l- z„)log 2 (l - z n ). 



(11) 



For nonuniform systems, the value of E^ n 
the site position n. At an eigenstate |\& 7 ). 
site-averaged von Neumann entropy 
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depends on 
we define a 



(12) 



were N is the system size. From the definition (|12[) . it 

shows that for an extended state that ib2 = -h? f° r all 

n vJV 

n, En v = -i log 2 jf - (1 - i) log 2 (l -jr)ajr log 2 N at 
A — > oo, and for a localized state that ip^ = 5 nn ° ( n° is 
a given site ) , = 0. In the paper all the values of E% 
and E1 n are scaled by i log 2 N. From the two examples, 
we know the scaled E^ is near 1 when eigenstates are 
extended, and near zero when eigenstates are localized. 
Henceforth, we omit "scaled" for simplicity. 

In order to analyze the influence of system parameters 
on the von Neumann entropy for all the eigenstates, we 
define a spectrum-averaged von Neumann entropy as a 
further gross measure, i.e., 
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(13) 



FIG. 2: The fidelity F(X, u; A , u ) and its contour map 
as functions of (A,u) at (a), (b) and (c) for 
(A ,u )=(l. 0,0.5), (2.0,1.5) and (3.0, 0.75), which 
corresponding to the system in the metal phase I, the 
metal phase III, and the insulator phase II, respectively. 



where M is the number of all the eigenstates. 



III. NUMERICAL RESULTS 



In numerical calculations, without loss of generality, 
we set k y — 0. As a typical case, (f> = (V5 — l)/2. In 
fact as is customary in the context of quasiperiodic sys- 
tem, the value of (f> may be approximated by the ratio of 
successive Fibonacci numbers — F m = F m —2 +-F m _i with 
Fq = F\ = 1. In this way, choosing <p — F m —i/F m and 
system size N = F m , we can obtain the periodic approx- 
imant for the quasiperiodic potential. We directly diago- 
nalize the eigenvalue Eq.([2]) at different values (A,/z) and 



get all the eigenvalues E 1 and the corresponding eigen- 
states |^ 7 ). From the formulas fTB"]) . we can obtain 
the fidelity F(X, u; Ao, uo), the FS \f, the site-averaged 
von Neumann entropy E^ and the spectrum-averaged 
von Neumann entropy (E v ), respectively. Henceforth, 
for simplicity we denote F to F(X, fx; Ao, fXo) ■ In all the 
figures the system sizes N is chosen to Fibonacci number 
987 unless specially stated. 



A. Fidelity 

In the metallic phase I, metallic phase III, and insulat- 
ing phase II, we choose (Ao, /zq)=(1.0, 0.5), (2.0, 1.5) and 
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(3.0,0.75) as examples, respectively. The corresponding 
fidelity F varying with parameters A and fj, are shown in 
Fig|2] At the same time, the contour maps of the fidelity 
are also shown. It shows that when parameters are at the 
same phase, the fidelity is near one; otherwise, the fidelity 
is very small. It is interesting that, though the phase 
I and III are both metallic phases and the correspond- 
ing wave functions are all extended, the fidelity is small 
when parameters are in the two phases respectively. This 
can be understood from the corresponding "classical or- 
bits" Hamiltoniar*2&: for Phase I, the contour lines of the 
Hamiltonian are extended in x direction but localization 
in the y direction, while for Phase III, the contour lines 
are extended in x + y direction but localized in the x — y 
direction. Therefore, the two phases are different. At 
the same time, these contour maps divide the parameter 
space to different regions, which is good agreement with 
the phase diagram shown in Fig[T] Comparing with the 
fidelity shown in Fig[5Ja) and (b), the fidelity in FigJ^tc) 
changes drastically with model parameters. It is because 
the band edge states in the insulating phase II may be 
localized in different regions of space and the overlap of 
these states may be large or small. 

At the three phase boundaries, we choose 
(Ao,Mo)=(l-0 ) 1.0), (2.0,0.5) and (3.0,1.5) as examples, 
which corresponding to the system at the boundaries 
between Phase I and III, Phase I and II and Phase III 
and II, respectively. The fidelity F varying with A and 
fi are shown in FigJ3] It shows that when parameters 
(A, ^i) and (Ao, /io) are at a same critical line, the fidelity 
is near one; otherwise, the fidelity is relatively small. It 
is interesting that if a point (Ao,/io) in the critical line 
between Phase I and II(Phase I and II, Phase II and 
III), the fidelity in the both phases is relatively large. 
Similar to that shown in FigJ^l these contour maps of 
fidelity also divide the parameter space to three regions, 
which is same as the phase diagram shown in Fig[TJ 

At the bicritical point (Ao, /j.o)=(2.0, 1.0), the fidelity 
F and its contour map as functions of A and /i are plot- 
ted in Figgl It shows that when (A, /i) = (2.0, 1.0), F is 
maximal and equal to one, when (A, fi) for the three crit- 
ical lines, F becomes relatively small, and when (A, fi) 
in Phase I, II and III, F becomes relatively smaller. All 
these certify that the bicritical point itself is different 
from others points, which is agreement with the conclu- 
sion that the bicritical point is a particular critical point 
as investigated in Ref^. At the same time, the contour 
of fidelity can reflect the phase diagram shown in Fig[TJ 



B. Fidelity Susceptibility 

According to the definition of FS in Eq. © , its values 
depend on g a b and a specific direction of parameter path 
q = q(\,/j,)H. The tangent unit vector (n , n* 4 ) which 
defines the direction may be different, though g a b does 
not depend on parameter paths. In the following, the 
FS for (n\n^) = (0,1), (1,0), and (1/VE, -2/VE) are 



1.0 




FIG. 3: The fidelity F(A,/i; Xq^o) and its contour map 
as functions of (X,/i) at (a), (b) and (c) for 
(A ,/i ) = (1.0,1.0), (2.0,0.5) and (3.0, 1.5), which 
corresponding to the system at the phase boundaries 
between I and III, I and II, and III and II, respectively. 



1.0- 




FIG. 4: The fidelity F(A,/i; \o,Ho) an d its contour map 
as functions of (A,/z) for (Ao, /xo)=(2.0, 1.0), which 
corresponding to the system at the bicritical point. 
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(c) 

FIG. 5: The logarithmic plots of fidelity susceptibility 
Xf and its contour maps as functions of (A, (i) from 
different parameter paths q = q(X, fi). The tangent 
units vector (n A ,n AI ) of paths equals to (a)(0, 1),( 
b)(l,0), (c)(l/V5,-2/%/5), respectively. 



shown in Figj5|a), (b) and (c), respectively, which are 
corresponding to that only \x , only A changes and both 
change simultanely. 

Fig[5] shows the Xf(\ m) are different when choosing 
different parameter paths. For Figf5]Ja), only the driv- 
ing Hamiltonian effects the values of xf- From the 
corresponding contour maps, the boundaries between the 
metallic phase III and other two phases are identified, i.e., 
there are sharp changes in xf at these phase boundaries, 
while for Fig[5f b), only H\ effects xf and the boundaries 
between the insulating phase II and other two phases are 
identified. The combination of the two contour maps in 
FigGUa) and (b) is consistence with the phase diagram 




FIG. 6: log (AE) varying with (A,/i), here AE is the 
gap between E\ and Eq at (A,/i). 



that shown in FigJTJ For FigO^c), both H\ and ef- 
fect xf and the corresponding contour maps itself can 
reflect the phase diagram. In FigGJa), (b) and (c), the 
varying of xf in the insulating phase II is not smooth, 
which is due to the gap between Eq and E 1 may be close 
to zero at some parameters(see Eqj9j). For this, the log- 
arithmic plots of the gap AE for the first excited state 
eigenenergy Ei and ground state eigenenergy E varying 
with (A, /i) are shown in Fig[Hl One sees that in Phase II, 
all the values of AE are very small and some almost are 
equal to zeros, therefore the fidelity F changes sharply at 
these parameter points. It is interesting that the contour 
maps of AE divide the parameter space to three regions, 
which is also consistence with the phase diagram shown 
in FigEEJ 

In order to study the critical behavior around criti- 
cal points(A c , /i c ), we study the finite scaling analysis of 
F S 13 ' 15 i 16 ' 17 and obtain the corresponding critical expo- 
nents. It has been found that these critical exponents and 
different scaling behaviors of FS can characterize the uni- 
versality classes of phase transitions^. Firstly, we study 
the transition between the metallic phase I and metallic 
phase III and choose the critical point (A c = 1.0, [i c = 
1.0) as an example. Near the critical parameter with 
the tangent units vector (n x = 0,n M = 1) of parameter 
paths, the xf is calculated for various system sizes N, 
which corresponding to the case that shown in Figj6ja). 
Along the parameter path, the FS reaches its maximum 
value XFmax at a certain position jJb max . The scaling be- 
haviors of XFmax and [iFmax - Mc are given in Fig^a) 
and (b), respectively, which shows that XFmax a N a and 
HFmax — . All j3 < 0, which means /iFmax tends 

to the critical point fi c in the thermodynamic limit. For 
the system sizes TV are chosen to the Fibonacci number 
F m with m = 31 + 1 and m 7^ 31 + 1 for integer I, the 
system sizes are divided to two cases 3 ^. It is found that 
a w 2.0 and w -2.0 for m ^ 31 + 1, while a = 4.9371 
and /3 = -1.5022 for to = 31 + 1. In FiglgJ the corre- 
sponding scaling functions are plotted. It shows that the 
exponent v w 1.0 for to 7^ 31 + 1, while v = 2.4718 for 
to = 31 + 1. Although a(/3)are different for m = 3/ + 1 
and to 7^ 31 + 1, the scaling relation a/v sa 2 is universal. 
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FIG. 7: The scaling behaviors of (a) XFmax and (b) 
HFmax — Mo respectively. The system sizes 
F 3l = 55, 233, 987(o),F 3i+ i = 89, 377, 1597(D) and 
F 3 i +2 = 144,610, 2584(A), respectively. At the same 
time, the corresponding fitted lines are also shown, 
respectively. 



The scaling relation is same as that for a one-dimensional 
asymmetric Hubbard model studied by Gu et al^. 

To understand the different behaviors between the sys- 
tems with N — F 3i+ i and N ^ F 3t+ i, we analyse care- 
fully the structure of the system. According to Fibonacci 
numbers F m = F m _ 2 + F m -i with F a = Fx = 1, F 3t 
and F 3 i + i are odd, which can be written as 2ki + 1 and 
2&2 + 1 with integers k\ and k 2 , respectively. For N — 



F 



3i+li 



F31 



2fci+l 
2k 2 + l 



Eq.©,-[ 
at the site n 



1 -I- i^W) 



and fi = 1, the hopping term of 

■ 2k 1 + l 2n+l 

-[1 + e 2 "^ 2 +n-] = o 
a bond between k 2 and (k 2 + l)th 



= k 2 , i.e 

site breaks. The system is divided to two segments. For 
N 7^ Fai+i, it docs not happen. This induces differ- 
ences between the energy spectrum of N = F 3 i + \ and 

a /•:<„■ . A 

Secondly, we study the transition between the metallic 
phase I and insulating phase II and choose the critical 
parameter (A c = 2.0, fi c = 0.5) as an example. Near the 
critical parameter with the tangent units vector (n A = 
1,71^ = 0), the FS xf is calculated for various system 
sizes N, which corresponding the case shown in Fig(5Jb). 
From Figs. l9"land[TT)l it is found that for all system sizes, 
55, 89, 2584, the exponents a, f3 and v are same and 
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FIG. 8: The finite size scaling analysis is performed for 
{XFmax ~Xf)/xf as a function N v {n - Umax)- 



the the scaling relation a/v « 2 is also obtained. We 
have studied the transition between the metallic phase 
III and insulating phase II, the results are similar and 
the relation a/v w 2 is also tenable. 



C. von Neumann entropy 

The von Neumann entropy has been found to be a 
suitable quantity to characterize the localization proper- 
ties of electronic state s 25 i 28 ' 29 . Fig fTTT a) and (b) show 
the site-averaged von Neumann entropy E% for the low- 
est edge states and the spectrum-averaged von Neumann 
entropy {E v ), respectively. The varyings of the two quan- 
tities with parameters(A, (i) are similar. E^((E V )) is near 
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FIG. 9: The scaling behaviors of (a) XFmax and (b) 
^Fmax — A c , respectively. The system sizes are 
55, 89, ...,2584. 




FIG. 11: The site-averaged von Neumann entropy E^ 

for the lowest edge states (a) and the 
spectrum-averaged von Neumann entropy (E v ) (b)as as 
functions of (A,/i), respectively. 
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FIG. 10: The finite size scaling analysis is performed for 
{XFmax - Xf)/xf as a function N v {\ - \ max ) for the 
system sizes are 55, 89, 2584. 



1 in the metallic phase I and III and relatively small 
in the insulating phase II. There are sharp decreases in 
E2((E V )) at phase boundaries. The contour maps of 
them divide the parameter space to three parts, which 
is consistence with the phase diagram shown in FigfTJ 

Conventionally, the inverse participation ratio (IPR) 
is often used as a measure of the wave-functions lo- 
calization length 39 . The larger is the IPR, the more 
delocalized the eigenstate is. It has found that the 



site-averaged von Neumann entropy E% increases ex- 
ponentially with the IPR 29 , i.e., E% can reflect the 
localization properties of electronic states. FigfTJ^a)- 
(d) show E~l varying with eigenenergy E 7 for (A, /i) = 
(1.0, 0.5), (3.0, 0.75), (1.0, 1.0) and (2.0, 0.5), which corre- 
sponding to the metallic phase, the insulating phase, the 
MMT and the MIT, respectively. In FigfTJ a), almost all 
the El are near 1, which means these states are delocal- 
ized. Comparing FigfT2Tb) with Fig[l"27a). all the El are 
small, which means that all eigenstates are localized. In 
FigfT2Tc) and (d), there coexist large, middle, and small 
E%, which means the eigenstates are critical. Though all 
eigenstates for the three phase boundaries and the bicriti- 
cal point are critical, the values of the spectrum-averaged 
von Neumann entropy (E v ) are different. Comparing to 
each other, the (E v ) for boundaries between the metal- 
lic phase I and III are large, for the bicritical point are 
middle, and for boundaries between metallic phase and 
insulator are small, which can be seen from Fig lllf b). All 
these indicate that, judging from the varying of von Neu- 
mann entropy with parameter (A,/i), the phase diagram 
can be completely characterized. 
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IV. CONCLUSIONS AND DISCUSSIONS 
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FIG. 12: The site-averaged von Neumann entropy 

varying with eigenenergy E 1 at 
(A, n) = (1.0, 0.5), (3.0, 0.75), (1.0, 1.0) and (2.0, 1.0) for 
(a), (b), (c) and (d), respectively. 



For the extended Harper model introduced in Ref. 36 , 
we have studied the fidelity between two lowest band edge 
states corresponding to different model parameters, the 
FS and the von Neumann entropy of the lowest band 
edge states, and the spectrum-averaged von Neumann 
entropy. All the three quantities can well characterize 
the rich phase diagram of the interesting model. 

In detail, firstly, the fidelity varying with parameters 
(A, n) for seven groups of fixing values (An, /xo) is studied, 
which corresponding to different phases, different phase 
boundaries and the bicritical point. When parameters 
are in the same phase or same boundary, the fidelity is 
near one, otherwise, it is small. There are drastic changes 
in fidelity when one parameter corresponding to phase 
boundaries. At the same time, the contour maps of fi- 
delity divide the parameter space to three regions, which 
is a good agreement with the phase diagram of the model. 
In fact, these conclusions are valid for arbitrary fixing 
values ( An, /in) in the parameter space. It indicates that 
the fidelity can well reflect the different phases and reveal 
different phase transitions. 

Secondly, the FS is studied and the finite scaling anal- 
ysis is performed for the MMT and the MIT. The con- 
tour maps of FS can well reflect the phase diagram. At 
the MMT, the critical exponents a(fi : v) for system sizes 
F m= 2.i+i and F m ^i + i are different, but the relation that 
vja w 2 is universal. At the MIT, the critical expo- 
nents for all system sizes are same and the relation that 
v/a w 2 is also tenable. 

At last, the von Neumann entropy is studied. It is near 
one in the metallic phase, while small in the insulating 
phase. There are sharp changes at phase boundaries. 
There are difference in the values of spectrum-average 
von Neumann entropy for the three phase boundaries 
and the bicritical point. The contour maps of von Neu- 
mann entropy is consistence with the phase diagram. All 
these indicate that the different phases and phase tran- 
sitions can be completely distinguished by von Neumann 
entropy 
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